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Mesons consisting of a light and a heavy quark are an 
interesting laboratory for the study of quantum chromo- 
dynamics, the theory of strong interactions. The role of 
the heavy-light quark mesons is similar to that played by 
the hydrogen atom in understanding quantum electrody- 
namics. The bound states of a & quark with either a light 
u or d quark are referred to as B mesons. The states with 
zero internal orbital angular momentum (L — 0) and spin 
parity J p = 0~ (B) and l - (B*) are well established pj , 
but the spectroscopy of the orbitally excited B states has 
not been well studied. For L = 1, the total angular mo- 
mentum of the light quark is j = | or j = |. With the 
addition of spin of the heavy quark, two doublets of states 
are expected: states with j — |, named Bq ( J = 0) and 
B[ (J — 1), and states with j — |, named Bi (J = 1) 
and B\ (J — 2). These four states are collectively re- 
ferred to as B**. 

Heavy quark effective theory Q predicts that the mass 
splitting within each doublet of a heavy-light quark me- 
son is inversely proportional to the heavy quark mass 
@, 0, II 0, H 0, II- The 3 = \ states are expected to 
decay to B^tt via an S-wave transition and to exhibit 
resonance widths in the range 100-200 MeV/c 2 [§]. The 
j = I states are expected to decay to B^ir via a D-wave 
transition and to have widths of 10 - 20 MeV/c 2 0, Q. 
This Letter focuses on the Bi and B% observed in Bit 
final states. The decay B\ — > Bit is forbidden by con- 
servation of angular momentum and parity, while both 



B* 2 



Bit and B^ — > B*n decays are allowed. Decays 



to a B* are followed by B* — > Bj, where the photon 
is not reconstructed in CDF due to its low energy. Be- 
cause of the missing photon, the measured Btt mass in 



Bi 



B*ir — > Bnj and B% — > B*ir — > Bnj events is 



lower than the B*ir mass by 45.78 ± 0.35 MeV/c 2 Q, 
resulting in an expected signal structure of three narrow 
Btt peaks for the J5i and £?|. 

Previous measurements of properties of the j ; = | B® 
and B*° mesons using inclusive or partially reconstructed 



decays did not separate the narrow states T3, U| or were 
limited by low sample statistics [HI]. Recentl y th e DO 
Collaboration resolved the B® and B^ masses hj. The 
superb mass resolution of the CDF II detector allows bet- 
ter precision and enables us to measure the B^ width. 
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Here, we present measurements of the masses of the B® 
and B^ states and the width of the B^ state. We recon- 
struct B**° in B + n~ and B* + ir~ decays, where the B + 
candidates decay into J/ipK + , D°ir + , and Z) 7r + 7r + 7r~ 
final states with J/ip — > H + /J.~ and D° — > K + ir~. 
Throughout this paper, any reference to a specific charge 
state implies the charge conjugate state as well. 

We use a data sample of events produced in pp colli- 
sions at ^/s = 1.96 TeV recorded by the CDF II detector 
at the Tevatron, corresponding to an integrated lumi- 
nosity of 1.7 fb . The components and performance 
parameters of CDF II 14j most relevant for this analysis 
are the tracking, the muon detectors, and the trigger on 
displaced vertices. The tracking system lies in a uniform 
axial magnetic field of 1.4 T. The inner tracking volume 
is instrumented with a layer of single-sided silicon mi- 
crostrip detectors mounted directly on the beam pipe at 
a radius of 1.5 cm, and 7 layers of double-sided silicon 
that extend out to a radius of 28 cm [15| . This system 
provides excellent resolution of the impact parameter, do, 
defined as the distance of closest approach of the track 
to the interaction point in the transverse plane. The 
outer tracking volume contains an open-cell drift cham- 
ber (COT) up to a radius of 137 cm [lf|. Muons are 
detected in planes of drift tubes and scintillators [T3| lo- 
cated outside the hadronic and electromagnetic calorime- 
ters. The muon detectors used in this study cover the 
pseudorapidity range \r]\ < 1.0, where 77 = — lntan(#/2) 
and 9 is the polar angle measured from the proton beam. 

A three-level trigger system selects events in real time. 
A dimuon trigger 14| requires two tracks of opposite 
charge that match track segments in the muon cham- 
bers and have a combined dimuon mass consistent with 
the J/ip mass. An extremely fast tracker at level 1 
(XFT) [11] groups COT hits into tracks in the trans- 
verse plane. A silicon vertex trigger at level 2 (SVT) [l9|] 
adds silicon hits to tracks found by the XFT, thus pro- 
viding better-defined tracks and allowing candidate selec- 
tion based on the impact parameter. A displaced vertex 
trigger [2(| requires two tracks each with a scalar trans- 
verse momentum, pt, greater than 2 GeV/c and with 
0.12 < do < 1 mm. Additionally, the intersection point 
of the track pair must be transversely displaced from the 
pp interaction point by at least 0.2 mm, and the pair 
must have a scalar sum Pt(1) +Pt(2) > 5.5 GeV/c. 

Decays B + — > J/tpK + are reconstructed from the 
dimuon trigger data while decays B + — * D° n + (n + n~ ) 
are reconstructed from the displaced vertex trigger data. 
In each decay, the tracks are constrained in a three- 
dimensional kinematic fit to the appropriate B + vertex 
topology with the J/ tp and D° masses constrained to the 
world average values [l[. Each track compatible with 
originating from the same interaction point as the B + 
and not used to reconstruct the B + is considered as a 
pion candidate, and its 4-momentum is combined with 
that of the B + candidate to form a B**° candidate. We 



5 



search for narrow resonances in the mass difference distri- 
bution of Q = m(B + n^)—m(B + )~rn 7T , where m(B + 7r~) 
and m{B + ) are the reconstructed invariant masses of the 
B + 7r~ pair and the B + candidate, and is the pion 
mass. 

The B + candidates are selected using independent ar- 
tificial neural networks for each of the three B + decay 
modes. The neural networks are based on the Neu- 
roBayes package [21(. For the decays B + — ► J/ipK + 
and B + — > D°tt + , we use the training and selection 
methods developed in Ref. [22j]. For the decay B + — > 



D ^ 



we closely follow the construction of the 



neural networks for the other two decays. To train 
this last neural network we use data from the region 
5325 < m(B + ) < 5395 MeV/c 2 as the background sam- 
ple and simulated B + events as the signal sample [23j |. 
The most discriminating inputs to the neural networks 
are pt{B + ), do{B + ), do of the kaon or pion with respect 
to the B + decay vertex, and the projected distance of the 
B + decay vertex from the primary vertex along the B + 
transverse momentum. We select approximately 51 500 
B events in the J/tpK + decay channel, 40 100 in the 
D°ir + channel, and 11 000 in the Z)°7r + 7r + 7r~ channel. 

To select B** a mesons, three additional neural net- 
works are trained on a combination of a simulated signal 
sample and real data for a background sample. The data 
for the background sample are taken from the entire Q 
range of to 1000 MeV/c 2 , which includes only a small 
contribution from the signal in the data. To avoid biasing 
the network training, the simulated events are generated 
with the same Q distribution as the data. The B**° 
neural networks use the same inputs as the B + neural 
networks, together with the kinematic and particle iden- 
tification quantities for the pion from the B**° decay. 
The most important discriminants are the px and do of 
the pion from the B**° decay vertex and the output of 
the B + neural network. 

For each B + decay channel, we require fewer than six 
B**° candidates in an event in order to enhance the 
signal-to-background ratio. The observed B**° signals 
are consistent for all three B + decay channels. Therefore 
we combine the B**° events for all decay channels and 
use this combined Q distribution to measure the B**° 
properties. We count the number of Monte Carlo signal 
events, Nmc-i an d the number of signal and background 
events in the data, N^ata, in the Q signal region of 200 to 
400 MeV/c 2 for a given cut on each of the three network 
outputs. We then optimize the B**° selection for each 
B + decay channel to maximize the combined significance, 
Nmc / V Ndata- The resulting combined Q distribution is 
shown in Figure [TJ 

The B**° signal structure is interpreted as resulting 
from the three signal processes B® — > B* + ir~, _B|° — > 
B+tt", and B* 2 ° -> B*+ir-, with B*+ -> B+7. The Q 
distribution for each signal process is modeled by a non- 
relativistic fixed- width Breit-Wigner function convoluted 




Q = m(B7u) -m(B)- m,, (MeV/c ) 



FIG. 1: 

m(B + 7r 



Distribution of the mass difference Q = 
~)—m(B + )—m n for exclusive B + decays. 
Curves are shown separately for the background, 
the B**° -> B M K reflections, and the three 
decays. 



with the detector resolution model. The resolution on 
Q is determined from simulation and modeled as a sum 
of two Gaussian distributions, a dominant narrow core 
and a broad tail with Q-dependent standard deviations 
of about 2 MeV/c 2 and 4 MeV/c 2 , respectively. The 
fraction of events in the broad tail is fixed to be 0.2. 

We perform an unbinned maximum-likelihood fit to 
the combined Q distribution, from which we extract the 
Q value of the B^ — > B + ir~ decay, the mass difference 
between the _B° an d B%° states, the width of the B% , 
and the number of events in each signal process. The 
following parameters in the fit are constrained to their 
values from either previous measurements or theoreti- 
cal predictions: the energy of the B* + decay photon, 
E{i) = 45.78 ± 0.35 MeV/c 2 fl; the ratio of the and 

Bf widths, ^§U- } = 0.9 ± 0.2 0); and the ratio of the 
i?2° branching fractions, w^rSio ~^ 



j 2 uLauuiwgu^uuiffl, BR(Bf->B' + 7r-) ~ ^±0.3 [ll|, 

consistent with the value measured in Ref. fl3j |. 

The background is modeled by a sum of two compo- 
nents, each being the product of a power law and an 
exponential function. We also expect reflections from 
B**° — ► B + K~ decays when the kaon is mistakenly as- 
signed the pion mass. The shape of the reflection in 
the Q distribution is determined in simulations of B**° 
states 22] and fixed in the fit. The normalization of the 



B**° is obtained by correcting the observed yield from 
Ref. 22j by a ratio of efficiencies to reconstruct a B**° 
decay as a B**° and B**°. In the B**° data sample we 
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expect 24 ± 12 B° ± events and 62 ±31 B*% events. These 
normalizations enter the fit as Gaussian constraints. 

Sources of systematic uncertainty on the mass differ- 
ence and width measurements include mass scale, mass- 
dependent signal efficiency, fit model bias, assumptions 
entered as Gaussian constraints in the fit, choice of back- 
ground and resolution models, and location and amount 
of B**° broad states. The systematic uncertainties are 
summarized in Table [U 

To determine the mass scale uncertainty, we recon- 
struct ip(2S) — > J/ijjiT + iT~ with J/ip — > which has 
a similar Q value as the B**° decays. We compare the 
measured Q to the world average [l| and take the differ- 
ence as the mass scale uncertainty. To evaluate the effect 
of signal efficiency with changing Q, we generate a large 
number of samples of the same size as the data, called 
pseudoexperiments, with the Q-dependent efficiency ob- 
tained from simulation. We then apply the default fit to 
the pseudoexperiments. 

Tests of the fit model on pseudoexperiments show a 
small fit bias on the B**° signal parameters, which is 
included as a systematic uncertainty. Signal parameters 
entered as Gaussian constraints in the fit contribute to 
the fit uncertainty. To determine their systematic contri- 
bution, we refit the data with these constrained parame- 
ters fixed. This fit returns the statistical fit uncertainties, 
which are subtracted in quadrature from the total fit un- 
certainties to obtain the systematic contribution. 

To estimate the uncertainties due to the choice of 
background and resolution models, we generate pseudo- 
experiments with varied background parameterizations 
or worse mass resolution. The background is also well- 
modeled by the sum of a broad Breit-Wigner function 
with the product of a power law and an exponential func- 
tion. From comparisons of the detector resolution in data 
and Monte Carlo for the ip(2S) sample, we expect the 
Monte Carlo to underestimate the resolution by no more 
than 20%. These pseudoexperiments are fit with the de- 
fault fit and the generating model. The distribution of 
the differences between these fit results is modeled by 
a Gaussian, whose mean is assigned as the systematic 
uncertainty. 

Possible effects of Bq and B[ decays on our background 
model are studied by adding two Breit-Wigner functions 
of identical width varied over the range 100—200 MeV/c 2 . 
The Q values of the states are independently varied in 
the range 240 to 360 MeV/c 2 , the region around the nar- 
row B**° peaks. We refit the data for various masses and 
widths of the broad states, with the normalizations of the 
broad Breit-Wigner functions as additional free parame- 
ters in the fit model. We then take the largest variation 
in the narrow B**° parameters from any configuration 
of broad states as the systematic uncertainty due to the 
B**° broad states. 

The result of the likelihood fit to the data is shown in 
Figure [TJ and we measure the following: 



TABLE I: Systematic uncertainties on the B**° parameter 
measurements. Each row corresponds to one source of sys- 
tematic uncertainty. The columns show the uncertainties for 
each of the three B**° signal parameters. Uncertainties are 
in units of MeV/c 2 . 



Source 


Q(Bf) 


r(B 2 *°) 


m(Bf) - m(B? ) 


Mass scale 


±0.2 




< 0.01 


Efficiency 


+0 

-0.03 


+0.4 
-0 


+0 
-0.3 


Fit bias 


+0 
-0.1 


+0.4 
-0 


+0 
-0.2 


Fit constraints 


+0.4 
-0.3 


+2.1 
-1.5 


+0.7 

-0.9 


Background 


+0.2 
-0 


+0 
-1.6 


+0.2 
-0 


Resolution 


< 0.01 


+0 
-0.4 


< 0.01 


Broad states 


+0.7 
-0.5 


+2.3 
-9.9 


+0.9 
-1.0 


Total 


+0.9 
-0.7 


+3.2 
-10.2 


+ 1.2 

-1.4 



m(B*°)-m(B+)-m w = 321.5t^(stat.) ±j^(syst.) 
MeV/c 2 ; 

m{Bf) ~ m(Bf) = M.^^stat.^-^syst.) 
MeV/c 2 ; and 

T(Bf) = 22.7±^(stat.)l5 2 2 (syst.) MeV/c 2 . 

The signal is consistent with theoretical predictions @,@|, 
and Gaussian-constrained parameters remain close to 
their input values, the largest departure being 0.4 stan- 
dard deviations. The numbers of events are N(B®) — 
503±£f, N(B*° -> B+n-) = 385tg, and N{B*° -> 
B* + tt~) — 351^45, where uncertainties are statistical 
only. Using the mass of the B + [l[ and the correlations 
between the fit parameters, the masses of the B® and B|° 
are m(Bf) = 5740.2±^(stat.)l^(syst.) MeV/c 2 and 

m{Bl) = 5725.3l2:2( stat -)±i'.5(sy st -) MeV/c 2 . With the 
current statistics the data are also consistent with con- 
taining only the B® and B% — > B + ir~ peaks. 

In summary, using the three fully reconstructed decays 
B+ -> J/^K + , B+ -► D°ir+, and B+ -► 5 7r+7r+7r-, 
we observe two narrow B**° states in the decays B® — > 
B* + ir~ and B^ — > B^ + tt^ . This is the most precise 
measurement of the narrow B**° masses to date. We 
have also measured the B^ width for the first time. 
There is some discrepancy between these measurements 
and those reported by the DO collaboration [13j, the 
largest being close to a 3 a difference in the mass splitting 
of the two B**° states. 
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